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^ ■ ABSTRACT 
O 
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. We present Spitzer Space Telescope IRAC and MIPS observations of the galac- 

^ ! tic globular cluster M15 (NGC 7078), one of the most metal-poor clusters with 

a [Fe/H] = -2.4. Our Spitzer images reveal a population of dusty red giants 
near the cluster center, a previously detected planetary nebula (PN) designated 
, K648, and a possible detection of the intra-cluster medium (ICM) arising from 

\0 [ mass loss episodes from the evolved stellar population. Our analysis suggests 

^ I (9 ± 2) X 10~^ Mq of dust is present in the core of M15, and this material has 

^ I accumulated over a period of ~ 10^ yrs, a timescale ten times shorter than the 

\Q [ last galactic plane crossing event. We also present Spitzer IRS follow up observa- 
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tions of K648, including the detection of the [Ne II] A12.81 /im line, and discuss 
abundances derived from infrared fine structure lines. 



Subject headings: globular clusters: individual (M15) — ISM: evolution — stars: 
^ I AGB and post-AGB — stars: mass loss — planetary nebulae: individual (K648) 

> 

^ : 1. INTRODUCTION 

Mass loss in evolved stellar populations affects the chemical evolution of the interstellar 
medium (ISM) and mass loss from individual stars governs post main sequence evolution. 
The amount and duration of mass loss that occurs in giant stars remains one of the most 
uncertain parameters in stellar evolution theory and the effect of these processes on the 
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inferences derived from stellar population models can be significant. Given the wide use of 
these models (e.g., in inferring stellar masses of high-redshift galaxies), understanding the 
mass losing process is vital for a range of problems in astrophysics. Although dust constitutes 
a small fraction of the total mass lost, it is frequently used as a marker as it is optically thin 
and its thermal emission is readily detectable in a variety of environments. 

Globular clusters (GCs), believed to have formed during the assemblage of the Galaxy, 
are coeval samples of stars at common, well-determined distances with nearly uniform initial 
compositions. GCs enable study of the chemical enrichment of the interstellar medium 
arising from mass ejection during the post-main sequence evolution of stars. Red Giant 
stars, especially those ascending the Asymptotic Giant Branch, are expected to develop 
winds that inject processed material into the intra-cluster medium (ICM) during post-main 
sequence evolution. These winds contain gas and solid phase materials, the latter in the 
form of dust grains that condense from the metals. However, detection of thermal emission 
arising from intra-cluster medium dust has been elusive, suggesting that the ICM in GCs 
is 100 to 1000 times less massive than expected from current stellar evolution theory and 
observations of mass-losing stars in clusters and in the solar neighborhood. 

The circumstellar environments of stars in the late stages of evolution, when most 

mass loss occurs, arc most effectively detected and studied in the infrared (IR). IR surveys 
conducted by IRAS, 2MASS, and ISO revealed populations of dust-enshrouded asymptotic 
giant branch (AGB) and red supergiant (RSG) stars in the Galactic bulge (van Loon ct al. 
2003), the Large Magellanic Cloud (LMC) (Zijlstra ct al. 1996; van Loon et al. 1997; Trams 
et al. 1999), and in galactic Globular Clusters (Ramdani & Jorissen 2001; Origlia et al. 2002). 
GCs are expected to contain dust from episodes of mass loss in red giant branch (RGB) and 
AGB stars. The IR excess of their circumstellar dust emission is expected to peak between 
20 /im to 30 /im, and thus photometry at wavelengths larger than 20 /im is necessary to 
estimate accurate dust masses lost by such stars. The amount of dust present in the intra- 
cluster medium (ICM) will vary depending on the cluster escape velocity, the time since 
last crossing of the Galactic disk where the ICM can be stripped away by the ISM, and the 
number of mass-losing stars. In general the dust in the ICM of GCs is expected to be ~ 10^^ 
to 10~^ Mq for most galactic clusters. Globular clusters have reasonably homogeneous (in 
age and metallicity) , well-understood stellar populations, so observations of ICM dust are 
reasonably straightforward to interpret to yield mass loss rates and duty cycles. 

Previous attempts to detect the ICM in GCs suggest that the ICM density is well 
below that expected from predictions of the mass loss input from RGB and AGB stars, even 
considering the low metaUicity of GCs. The lowest S-cr upper limits to the ICM mass for 70 K 
dust are ~ 6 x 10~^ Mq (Hopwood et al. 1999). Detecting thermal emission from the elusive 
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ICM in GCs is observationally challenging. Origlia et al. (2002) reported ISO observations 
of the IR thermal emission from the winds of individual RGB stars in six massive GCs 
(47 Tuc, NGC 362, Omega Cen, NGC 6388, M15, and M54) showing that, in those systems, 
stellar winds from these stars are enriching the ICM. Though thermal emission from the ICM 
material might be expected to be detectable, many attempts to do so with IR and millimeter 
observatories have produced only a single secure detection of ICM dust, a tentative (3.5-(t) 
detection of thermal emission in the core of the metal-poor GC M15 (Evans et al. 2003). 
Overall this result suggests that the ICM dust in GCs is significantly less massive than 
expected from current stellar evolution theory and observations of mass losing stars in GCs 
and the solar neighborhood. 

The causes of the paucity of ICM emission have been proposed to be ram-pressure 
stripping of ICM gas during Galactic plane passage, blowout by nova explosions, fast winds 
from the stars themselves, radiative ejection by the sheer luminosity of cluster stars, and 
continuous ram-pressure from hot gas in the galactic halo. However, the dominant state of 
the ICM is unclear. Sensitive searches for neutral H in the ICM at radio wavelengths have 
yielded upper limits < 0.1 M0 (Birkinshaw et al. 1983), with a possible detection of ~ 200 
Mq in NGC 2808 (Faulkner et al. 1991) and a Sa detection of 0.3 Mq in M15 (van Loon et 
al. 2006). Perhaps much of the ICM is ionized, as suggested by the high electron densities 
measured from pulsar timing in 47 Tuc (Preire et al. 2001). However, Ha searches have as 
yet been unsuccessful. 

ICM dust grains in radiative thermal equilibrium should attain temperatures of 50 K 
to 80 K because of the high energy density of starlight within a GC (Forte et al. 2002) and 
thus be detectable as an IR excess (above the photospheric emission) in GCs at mid- to 
far-IR wavelengths. Here we present observations of the galactic GC M15 with the Spitzer 
Space Telescope (Werner et al. 2004), whose instrumental sensitivity enables detection of 
dust masses as low as 4 x 10~^ Mq (assuming a population of grains radiating in thermal 
equilibrium, with integration times down to the background confusion limit) and therefore 
permits an unequalled opportunity to search for and set stringent limits on the "missing" 
ICM. 



1.1. M15 

M15 (NGC 7078), with [Fe/H] = -2.4 (Sneden et al. 1997), is one of the most metal 
poor GCs. It is a well-studied cluster, as it is home to the first planetary nebula (PN) 
discovered in a GC (K648, also designated as Ps-1, Pease (1928); Howard et al. (1997); 
Alves et al. (2000)) and to the first GC Low-Mass X-ray Binary source (X2127-M19, Auriere 
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et al. (1984); Charles et al. (1986)). At least eight millisecond pulsars are also associated 
with the cluster (Kulkarni & Anderson 1996). M15 is generally believed to be a core-collapse 
GC, with a small, dense core containing approximately 4000 Mq (Phinney 1996). 

Properties of M15, reproduced from Hopwood et al. (1999), are listed in Table 1. Up- 
dated values as hsted by Evans et al. (2003) include the escape velocity Vesc from Webbink 
(1985), and the time Tc since the last plane crossing from Odenkirchen et al. (1997). The 
reddening and the distance are updated from Schlegel et al. (1998) and McNamara et al. 
(2004), respectively. Galactic coordinates for M15 are / = 65.01° and b = —27.31°, placing 
the cluster ~ 4.5 kpc south of the galactic plane. 

The total dust mass expected in a GC can be estimated using the following equation: 

lQ[Fe/H] 

Mdust = Nhb SM (1) 

Thb 100 

where thb is the Horizontal Branch (HB) lifetime, Nhb is the number of HB stars, [Fe/H] 
is the cluster metaUicity, and SM is the dust mass lost from each star at the tip of the RGB. 

The factor of 100 is the solar gas-to-dust ratio, which is scaled for the metallicity of M15 by 
adding the [Fc/H] factor. Using values typical of population II stars and this relationship, 
the expected dust mass in M15 has been estimated to be 3.7 x 10~^ Mq by Evans et al. 
(2003) and 2.0 x lO'^ Mq by Hopwood et al. (1999). 

Ml 5 is also home to the PN K648. K648 was the first globular cluster PN discovered, 
and subsequently it has been extensively studied at UV, optical and IR wavelengths to 
determine the chemical composition of the ejecta nebula and parameters of the central star 
(Barker et al. 1983; Adams et al. 1984; Howard et al. 1997; Bianchi et al. 2001; Dinerstein et 
al. 2003; Garnett & Lacy 1993). The study of post-AGB stellar evolution in old metal-poor, 
low-mass stars in GCs or the galactic halo can be greatly enhanced if the by-products of 
stellar nucleosynthesis can be measured. Enriched material produced in the RGB and AGB 
stages of stellar evolution are dispersed into outer layers of the stellar system and subsequent 
mass loss processes lead to the formation of a PN. Few PNe in the galactic halo population 
have been identified and only four of these PNe, including K648, BB-1, DdDM 1, and H4-1 
are associated with globular clusters (Jacoby et al. 1997). Below, we present findings derived 
from a 5-15 /xm IR spectrum of K648. 
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2. OBSERVATIONS 

Image observations of the GC M15 were obtained on 2004 October 29 UT with the 
Multiband Imaging Photometer for Spitzer (MIPS) (Rieke et al. 2004) camera through the 
24 /xm and 70 /xm filters and with the Infrared Array Camera (IRAC) (Fazio et al. 2004) 
at 3.6, 4.5, 5.8, and 8 /xm conducted as part of the Gehrz Guaranteed Time Observing 
Program (GGTOP, PID 124). Raw data were processed with the Spitzer Science Center 
(SSC) pipeline version Sll.4.0. To avoid saturation in IRAC, High Dynamic Range (HDR) 
mode was implemented, using 0.4 s, 10.4 s, and 98.6 s exposures. For each IRAC channel, 
55 frames at seven dither positions were obtained with each HDR exposure time. With the 
MIPS camera, 596 dithered frames were taken at 24 /im along with 256 dithered frames at 
70 iim. 

Spitzer IRS spectra of the PN and other red sources detected in the IRAC images 
have been obtained as part of a follow on program using GGTOP time and these data are 
discussed below. Table 2 summarizes the observations discussed in this paper. 

2.1. IRAC 

The IRAC Basic Cahbrated Data (BCD) images were post-processed to correct for var- 
ious instrumental artifacts and were mosaicked using the 2005 May 09 version of the SSC 
Legacy MOPEX software (Makovoz & Marleau 2005). The MOPEX cosmetic correction was 
used to eliminate mux-bleed and column pull-down as described by the IRAC Data Hand- 
book, v3.0.^ The background matching correction was used to minimize pixel differences 
in overlapping areas of the mosaics. The MOPEX mosaicker also eliminated cosmic rays 
and other outliers in the data. The final IRAC mosaics are not subsampled, thus the pixel 
size of each final tiled image is 1.22" per pixel covering an area ~ 50 square arcminutes 
around the core (RA[J2000] = 21^:29™:58".38; Dec[J2000] = +12°:10':00".6) of M15 plus 
an off-field region of the same size 7' to the north or south. A three color image combining 
3.6 /xm, 4.5 /xm and 8 //m is presented in Fig. la. The diffuse ICM is not detected in the 
IRAC mosaics, even at 8 ^m. However, the planetary nebula (K648) and several other dusty 
stellar objects become quite prominent at 8 /im. 

Point source photometry was conducted using several stellar extraction routines includ- 
ing GlimpsePhot^, DAOphot (Stetson 1987), and the astronomical point source extraction 
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(APEX) tool contained in the SSC MOPEX package. Severe stellar crowding towards the 
cluster core in the IRAC bands made reliable point-response-function (PRF) photometry 
challenging. The best photometric results (95% agreement to 2MASS K-band fluxes) were 
obtained by performing PRF fltting on the BCD images with APEX. Array location depen- 
dent photometric correction weights were applied before PRF fltting to minimize systematic 
error in the point source extraction flux. The images were then corrected with the back- 
ground matching routine, and outliers were eliminated by the mosaicker. 

PRFs were created for 3.6 fira and 4.5 fim data using the PRF estimate routine provided 
in APEX. At least twenty stars were chosen to make each PRF in areas with the smallest 
amount of crowding. The PRFs provided by the Spitzer Science Center (SSC)'^ that were 
made in-flight in January 2005 for 5.8 /im and 8 /xm are sufficiently over-sampled, so there 
was no need to create new PRFs for these channels. Probable point sources with fluxes at 
least 4-(7 above the background were then identifled and selected for extraction on coadded 
images that had been corrected for array distortion. Final flux extractions were performed 
at these world coordinates positions on the corrected BCD images using PRF subtraction. 

After applying the appropriate photometric color correction (as discussed in the IRAC 
Data Manual, v3.0) and comparing APEX 3.6 iim fluxes with 2MASS K-band fluxes (Skrut- 

skie et al. 2006), the average [K - 3.6 /im] color is that of a KO giant star. If we then assume 
that the majority of the stars we detect with Spitzer are KO giant stars, then our fluxes are 
consistent with 2MASS K-band fluxes to within approximately 5%. To further check our 
photometric results, we compared the fluxes from thirty of the stars we detect to K-band 
fluxes from Cohen et al. (2005). The median [K-3.6 /xm] color for these sources is 0.09, which 
is also consistent with KO giant stars (See Fig. 2). 

Color-magnitude diagrams (CMDs) derived from the IRAC point source photometry 
arc presented in Fig. 3. An RGB is clearly evident with the tip occurring near F3 g = 14 
mJy Stars with 3.6 /xm flux > 11.2 mJy are possibly saturated in the long exposure frames, 
so only medium exposure (10.4 second) BCDs were used for the PRF subtractions and flux 
extraction. Flux densities for stars with colors redward of the RGB are listed in Table 3 (See 
§3.1). 



^http://ssc. spitzer. caltech.edu/irac/psf.html 
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2.2. MIPS 

The MIPS Data Analysis Tool (DAT, Gordon et al. (2005)) version 2.96 was used to do 
the basic processing and final mosaicking of the individual MIPS images. In addition, extra 
processing steps on each image were carried out before mosaicking using programs written 
specifically to improve the removal of MIPS detector instrumental signatures. At 24 /im 
the extra steps included readout offset correction, scan mirror dependent flat flelds, a scan 
mirror independent flat fleld, array averaged background subtraction, and exclusion of the 
bias boost images. At 70 the extra steps were column average subtraction and pixel time 
filtering both with the exclusion of the regions around the bright sources. The pixel sizes 
of the final mosaics are 1.245" and 4.925" for 24 fim and 70 /xm, respectively. The 24 /xm 
mosaic covers a 90 square arcminute area and the 70 /xm mosaic covers a 5.0 arcminute x 
10.0 arcminute each centered on the core of M15. 

The MIPS mosaics (Fig. lb, c) show a possible ICM detection at 24 fim and 70 fim. 
The 24 fim image shows two high surface brightness patches of material, both offset from 
the core by ~ 17" towards the west (IRlb and IR2). The 70 fim image shows only one 
high surface brightness patch that is likely an unresolved image of both 24 /xm regions. 
The integrated fluxes of the ICM detections at both wavelengths were determined using 
basic aperture photometry in IDL v6.0. A 130.5" square aperture centered at RA[J2000] 
= 21^:29™:56".50; Dec[J2000] = +12°:09':53".32 yields fluxes of 159.4±0.1 mJy at 24 //m 
and 691.2±61.0 mJy at 70 /im. The 70 /im flux agrees well with that found by Evans et 
al. (2003) with ISO observations using the same aperture size and position (Table 4). Also 
visible in the MIPS mosaics are the planetary nebula at 24 fim and several dust enshrouded 
stars, some of which are also visible in the 8 /xm image, that may or may not be associated 
with the cluster. 

Point source photometry was performed on the 24 /xm mosaic using StarFinder (Diolaiti 
et al. 2000), which is well suited for the stable and well sampled MIPS 24 fim PSF. A 
STinyTim (Krist 2002)'^ model PSF with a temperature of 100 K and smoothed to account 
for pixel sampling was used for the stellar extractions. It has been shown that the smoothed 
STinyTim PSFs are exceUent matches to observed MIPS 24 /im PSFs (Engelbracht et al 
2006). CMDs comparing 24 //m to 3.6 iim and 8 /xm are presented in Fig. 4, and the 
uncertainties in the IRAC and MIPS photometry are summarized in Fig. 5. 
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2.3. IRS 

Associated with our imaging programs, observations of the PN K648 (Pease 1928; 
Howard et al. 1997; Alves et al. 2000) in M15 were obtained with the Spitzer Infrared 
Spectrograph (IRS) on 2005 November 17.77 UT using the short wavelength (5-15 /xm), low 
resolution module (SL) in staring mode. All observations utilized the IRS blue peak up array 
at the target position of the PN, RA[J2000] = 2l'^:29™:59".41; Dec[J2000] = +12°:10':25".70, 
and the entire Ha nebulosity of K648 (of. Fig. 2 of Alves et al. (2000)) was contained within 
the spectrograph slit. The slit dimensions are 57" x 3.6", and the slit was oriented 18.92 
degrees west of north. The SL spectroscopic astronomical observing templates consisted 
of 5 cycles of 60 second ramps. IRS BCDs were processed with version 13.0.1 of the IRS 
pipeline. A description of the IRS instrument and its operation is available in Houck et al. 
(2004). Details of the calibration and raw data processing are specified in the IRS Pipeline 
Description Document, vl.O.^ 

Post-pipeline processing was conducted to remove instrumental artifacts, perform back- 
ground subtractions and to combine extracted spectral segments. Fatally bad pixels were 
interpolated over in individual BCDs using bad pixel masks provided by the SSC. Multiple 
data collection events were obtained at two different positions on the slit using Spitzer 's 
nod functionality. The two-dimensional BCDs were differenced to remove the background 
flux contribution and then the data were extracted with the Spitzer IRS Custom Extractor 
(SPICE) (vl.l-betal5)^ using the default point source extraction widths. The extracted, 
background corrected data were combined using a weighted linear mean into a single out- 
put data file and clipped at the 3a level. At the time of reduction, the errors generated 
by the SSC pipeline were not reliable enough for sound interpretation and so the errors 
were estimated from the standard deviation of the flux at each wavelength bin. Where 
there were less than three points in a wavelength bin, the error is the quadrature sum 
of the errors in the files. The resultant spectrum is presented in Fig. 6 and derived hne 
fluxes are summarized in Table 5. The spectral lines were fitted using a least squares Gaus- 
sian routine that fits the line center, line amplitude, continuum amplitude and the slope 
of the continuum. The full-width half-maximum was fixed at the resolution limit of the 
low resolution module. Integrating the flux over the 8 fim IRAC bandpass yields a flux 
that agrees with the IRAC flux within the uncertainty limits. The strongest line in the 
mid-IR spectrum is the [Ne II] A = 12.81 fim line, followed by hydrogen recombination lines 
H I 6 — 5 (PfQ;)A = 7.46 /xm and H I 7 — 6 (Huq;)A = 12.37 /im. Emission from S^~^ is evident 



'"'http: / /ssc. spitzer. caltcch.cdu/irs/dh/PDD.pdf 
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in the spectrum, although the fit to the hne flux of the [S IV] A = 10.52 //m is of marginal 
signal-to-noise (S/N ^ 2.3), while no [Ar III] A = 8.99 /xm is seen. Our detections of the 
mid-IR neon and sulfur lines are the first reported in the literature for K648. Abundance 
estimates derived from these forbidden lines are discussed in §3.2. 

3. DISCUSSION 

Our Spitzer images of M15 for the flrst time clearly show both the stellar dust producers 
and the ICM dust (Fig. 1), allowing a direct comparison to be made between the dust 
injection and dust survival rates. The brightest source of 70 /xm emission, IRla, is blended 
at 70 /im but visible as separate objects in the MIPS 24 ^m map, IRlb and IR2. These 
sources arc completely invisible on the IRAC maps, even at 8 fim. A three color image 
of 8, 24, and 70 /^.m in which the 8 and 24 /xm images are degraded to match the 70 /xm 
resolution is presented in Fig. 7. This figure illustrates that IRla is not unresolved stellar 
emission, but is a starless dust cloud(s) that are likely to be of an intra-cluster nature. The 
next brightest objects at 70 /im are a pair of sources, IR3a and IR3b, that are situated at 
the fringes of the cluster and may not be physically associated with the cluster. IR3b was 
previously detected by 2MASS, but there are no known previous detections of IR3a. The 
probability of detecting non-member red sources in the field are sufficiently small to suggest 
that these sources are associated with the cluster. The only other 70 fim source, 1R4, is also 
located on the fringes of the cluster. Radial velocity measurements from Pilachowski et al. 
(2000) confirm its membership. Table 3 lists fiuxes for IR3a, IR3b, IR4, and other possibly 
dusty sources, described in §3.1. All three of these sources have mid-IR colors that are 
consistent with post-AGB stars (Fig. 4; Groenewegen (2006)). In addition to these sources, 
the planetary nebula K648 (Ps 1) is also detected in all IRAC bands and at 24 //m. The 
strong 24 /im dectection is likely due not only to dust, but also to line emission from [O IV] 
at 25.88 fim and/or [Fe II] at 25.98 fim. 

3.1. Mass-Losing Stars 

Figure 8 shows the locations of the mass-losing AGB stars in M15. The coordinates 
and IRAC fiuxes of these sources are listed in Table 3. These stars were identified by their 
locations on the IRAC CMDs (Fig. 3), redward of the RGB. We find 24 mass-losing/dust- 
enshroudcd stars and consider this a lower limit (due to potential source confusion) of the 
total number in the cluster. These stars also fall just redward of the RGB in the MIPS CMDs 
as well, which suggests that they could be post-AGB stars. Their IRAC colors, however. 
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indicate that it is more likely that they are AGB stars that are approaching the end of the 
AGB phase of their evolution (Groenewegen 2006). The stars blueward of the RGB can be 
explained by the absorption in the fundamental bands of CO in the 4.5 /xm band and SiO in 
the 8 fim band. The red, mass-losing stars in M15 populate an uneven spatial distribution 
about the core of the cluster as projected on the sky (Fig. 8). These stars are offset from the 
core in the same sense as IRla, IRlb, and 1R2, and the distribution is less cusped than the 
visual hght. Since M15 is 13.2 Gyr old (McNamara et al. 2004), most stars currently on the 
AGB have a zero-age MS mass of ~ 0.8 Mq. These stars will soon end their hves as white 
dwarfs of ~ 0.5 Mq, having lost approximately 0.3 Mq during their post-main sequence 
evolution. The loose spatial distribution of this population could therefore be due to mass 
segregation, in which lower-mass stars are displaced to the outer regions of the cluster due to 
their high velocities, leaving a preferential concentration of high-mass stars near the center 
of the potential well of the cluster (Spitzer 1987). 

The MIPS CMDs (Fig. 4) identify a population of approximately 23 post-AGB stars 
(Groenewegen 2006). The colors of lR3a, lR3b, and IR4 are similar to the colors of this 
population. These stars are bright at 24 fim, but they are not detected above the background 
at 70 fim (Table 3). They are distributed around the cluster center at an average radius of 
approximately 3.3', and their positions are biassed towards the southern side of the cluster, 
near IR3a, IR3b, and IR4. Mass segregation is most hkely responsible for their locations at 
the fringes of the cluster. 

Although M15 is one of the most massive galactic globular clusters, it is surprising to 
find many dusty objects, as the cluster metallicity of [Fe/H] = —2.4 places it amongst the 
most metal-poor galactic globular clusters. Dust forms from metal condensates and it is 
difficult to understand how dust grains can form at such low metal abundances. It is likely 
that the metals condensing to form dust are produced in the stars themselves and brought 
to the surface near the end of their evolution. The fact that dust production does not 
seem to be inhibited at metallicities ^1% solar implies that stellar mass loss must already 
have contributed dust very early on in the evolution of the Universe. Dust observed at 
high redshift is usually believed to originate in supcrnovac explosions that result from core 
collapse in massive stars, but our observations suggest that at least some of the dust formed 
within the first few 100 miUion years may have been produced by stars of only a few solar 
masses (e.g. MS hfetime from models by Vassihadis & Wood (1993)). 

Various investigators (Daulton et al. 2002; Amari et al. 2001) suggest that dust grains 
can form more easily at low metallicity in carbon stars, as these stars produce carbon them- 
selves. With less initial oxygen abundance to start with, it is easier for these stars to achieve 
C/0 > 1 in their carbon shells. With excess carbon available (after locking up equal amounts 
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of C and O in CO) , formation of carbon-chain molecules from which dust grain condensation 
can proceed. However, although there is evidence for high molecular abundances in stellar 
atmospheres at low mctallicity, dust production in these stars may be less efficient due to a 
lack of SiC seeds (Sloan ct al. 2006). Furthermore, van Loon (2006) suggests that the low 
optical depth in carbon stars in the Magellanic clouds points against large dust-to-gas ratios 
at low metallicity. 

For low-mass oxygen-rich stars, one would not expect this evolution to lead to effective 
dust production, as dredge- up only increases the C/0 ratio and does not facilitate the 
formation of oxygen-rich dust grains. In any case, for oxygen-rich stars it is well established 
that nucleation sites must be available to condense dust grains onto (Jeong et al. 1999). 
These seeds are hkely to be TiO or similar which include secondary elements that are not 
produced by the star itself. Other s-process seeds, such as Zr, may be dredged-up in the 
atmospheres of these stars, but the limiting factor for dust production with such seeds is the 
oxygen abundance, as more oxygen will be locked into CO after dredge-up. Unfortunately, 
our observations do not allow us to draw any conclusions as to the abundances of secondary 
elements in the mass-losing stars in M15. 

The suspicion is that stars do adjust their structure until they finally can shed their 
mantles as demonstrated by K648 in M15. From the analysis of spectra of GC giant stars, 
a picture is emerging in which metal-poor stars do become very cool while nonetheless 
exhibiting early-type spectra because the low metal abundances give rise to weak absorption. 
They may not form much dust, but they may still be able to form enough of it to drive a 
wind. The mid-IR spectrum of 47 Tucanae VI suggests typical silicate dust grains and a 
typical mass loss rate of 10~^ M© yr"^ (van Loon et al. 2006). 

3.2. IR Spectrum of K648 

Although UV and optical line ratios derived for K648 from previous investigations pro- 
vide constraints on relative abundance ratios, several important a-capture elements, such as 
S, Ar, and Ne, have ground configurations that produce only IR fine-structure. If lines from 
these ions are not observed and introduced into abundance models, the total abundance of 
these elements becomes uncertain (Dinerstein et al. 2003). Often, measuring emission fine 
flux from [Ne II]A12.81 //m, [S IV]A10.51 //m, and [Ar III]A8.99 //m is observationally chal- 
lenging from the ground. However, the sensitivity of the Spitzer IRS affords an opportunity 
to set stringent limits on the emission line fiux of these ions resulting in better constrained 
estimates on derived abundances. In addition, new radiative transition rates and collision 
strengths are now available as a result of the IRON project (Hummer et al. 1993) to improve 



- 12 - 



the accuracy of derived abundance ratios. Below, we discuss the analysis of our Spitzer 
measurements (Table 5) of the [S IV] and [Ne II] lines in K648 and present a reanalysis 
of the S/0 and Ne/0 ratios with contemporary atomic parameters using a simple model. 
Undertaking a full photoionization analysis of K648 (i.e., Howard et al. (1997)) using the 
Spitzer IR line fluxes is beyond the scope of this paper. 



3.2.1. H(3 

Accurate extinction correction of UV and optical lines with respect to H/3 are required 
for abundance determinations. Interstellar extinction estimates to K648 range from 0.2 < 
Ay < 0.6 (Garnett & Lacy 1993). However, at mid-IR wavelengths, extinction is minimal 
( ^ 0.02 mag), and the Pfa and Hua lines detections (Fig. 6) enable us to estimate the emitted 
H/3 flux directly. Assuming Case B (Osterbrock 1999), and adopting a Tg = 12,500 K 
and density of 10^ cm-^ (Garnett & Lacy 1993) we derived F{B.(3) = (1.91 ± 0.30) x 
10-^2 gj-gs ^j^2 g-i F(H/3) = (1.13 ± 0.36) x lO'^^ gj-gg ^^^2 g-i ^^ie Pfa and 

Hua lines respectively using the intrinsic hydrogen emissivity ratios tabulated in Hummer 
& Storey (1987). These two estimates are in reasonable agreement (within the formal error) 
with each other, considering there is an absolute photometric uncertainty of ~ 5% (Houck 
et al. 2004) between the two spectral orders (5.2 — 8.7 fim and 7.4 — 14.0 yum). Our average 
F(H/9) of (1.52 ±0.23) x 10^^^ ergs cm^ s~^ is in good agreement with previous observational 
estimates, especially those obtained with large apertures (Garnett & Lacy 1993), and we will 
adopt this value in our abundance analysis. 



3.2.2. Neon 

The ratio of [Ne II]A = 12.81 fim line flux to our derived average value of H/9 was used 
to estimate the Ne^/H^ abundance. We have adopted a collisional strength T(T) = 0.283 
(appropriate for Tg = 10^ K, Saraph & TuUy (1994)) and an Aif value of 8.59 x 10~^ s~^ 
from the NIST database and assumed an electron density of A^e = 1-7 x 10^ cm~^ (Garnett 
& Lacy 1993). Rate coefficients, qfi (Hummer et al. 1993) and the population levels were 
computed assuming Ne"^, a 2p^ ion, is a two- level atom. Following Rank et al. (1978), we 
define the relative abundance of Ne'^/H+ as 

Ne+ ^ {in jHp/N,N,)N, /([Ne II]) 

H+ huf.Afff I{H(5) ' ^ > 
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where the H/3 volume emissivity, ^TrjHp/NeNp, is 1.0301 x 10"^^ erg cm^ interpolated 
for the assumed Ng (Hummer & Storey 1987), // = 9.86 x 10~^ is the population ratio 
of the upper to lower state, and I([Ne II]), the observed neon line flux, is (1.74 ± 0.08) x 
10-13 erg cm-2 (TableS). This yields a ratio of Ne+/H+ = (1.53±0.21) x 10"^ The rela- 
tive abundance of Ne^+/H+ was estimated from optical measurements of the [Ne III]A3869 A 
and [Ne III]A3967 A lines (Adams et al. 1984) and de-reddened assuming c = 0.12 and a 
Seaton (1979) extinction law. Using the [Ne III] to H/3 ratios (see Table 6), the relative pop- 
ulations were computed using a multi-level atom program incorporating the best available 
Aij{s~^) and coUision strengths available from the NIST database and the literature (Moore, 
Hester, & Dufour 2004), which is similar to code originally described by Shaw Sz Dufour 
(1995). Summing over ah ions, we find a total Nc/H= (2.39 ± 0.27) x 10"^ and a Ne/0 
= 0.48 ± 0.14 adopting an 0/H = (5.0 ± 1.3) x 10"^ (Pena et al. 1993). 



3.2.3. Sulfur 

The relative sulfur abundance was computed in a similar manner as the neon abundance 
described in §3.2.2, although all relative populations were determined using the (Moore, 
Hester, & Dufour 2004) code. The Spitzer observation of [S IV] was used to determine the 
S^"'"/H+ population, while the optical fluxes reported by (Barker et al. 1983) corresponding to 
B.(3 (assuming j{Ha)/j{H(3) = 2.81; Osterbrock (1999)) were used to estimate the relative 
g2+^-[^+ g^^^ S+/H+ abundances (Table 6). We find an upper limit to the total S/H of 
(4.28 ± 1.28) X 10"®, which is ^ 2.5 lower than that inferred by Garnett & Lacy (1993). 
Adopting an 0/H = (5.0 ± 1.3) x 10"^ (Pena et al. 1993), yields a S/0 = (8.56 ± 3.39)"^. 



3.2.4. Abundance Comments 

Our new estimates of the [S/0]< —2.64 and [Ne/0] = +0.54 confirm that K648 is 
under-enriched in S as compared to O (Garnett & Lacy 1993), while Ne/0 is enhanced with 
respect to solar. The Ne enhancement is seen in other halo population PNe, such as BB-1 
(Dinerstein et al. 2003). Garnett & Lacy (1993) argue that contamination of He-burning 
products by a-captures at high temperature could account for enhanced Ne. Bianchi et 
al. (2001), based on analysis of Hubble Space Telescope (HST) Faint Object Spectrograph 
(FOS) spectra of the central star, suggest that the nebular shell was ejected by a low- 
mass He-burning progenitor that has subsequently undergone a late thermal pulse, perhaps 
similar to the evolution of objects akin to FG Sge (Gehrz ct al. 2005). Our derived neon 
abundance using the fine structure line, suggests that dredge-up from the stellar core may 
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be an important mechanism to pollute the expelled nebular materials from slowly evolving, 
young PNe. 



3.3. Dust in the ICM 

Assuming that the diffuse emission detected in the Spitzer images arises from ICM ma- 
terial, we can compute a dust mass using the observed SEDs. The approximate temperature 
of the ICM dust was derived by least-squares fitting a greybody model to the SED using 
fluxes measured in a 130.5 x 130.5 square arcsccond aperture centered at the same right as- 
cension and declination coordinates in images at all wave bands. Our choice of aperture size 
is equivalent to that used by Evans et al. (2003) to facilitate direct comparison. The large 
wavelength range of the Spitzer SED enables us to distinguish between the contribution of 
a stellar blackbody (peaking near 0.64 /i,m and dominated by KO stars), and that of thermal 
radiating dust that generates an IR excess at wavelengths greater than 24 /im. Use of a two- 
component model, incorporating a stellar blackbody that peaks near 4699 ± 58 K and a dust 
blackbody that peaks near = 70±2 K, gives a rough fit to the data (Fig. 9; Table 4). The 
fit produces a large reduced value of 5.26 suggesting that the integrated flux within our 
aperture sums the emission from stars of many disperate spectral types not well-represented 
by a simple, single emissivity and temperature blackbody. Another source of uncertainty in 
the fit is the effect of crowding in 2MASS data. This could lead to oversubtraction of the 
background and overestimation of the flux densities (van Loon et al. 2005). 

The mass of the ICM was determined using the methodology described in Evans et al. 
(2003). We assume that the dust is optically thin, which yields the following expression 

M 

where Dkpc is the distance to M15 in kiloparsecs (Table 1), Ki, is the dust absorption coefficient 
in cm^ g^^, B{v^Tfj) is the Planck function in cgs units, and is the dust temperature, 
at 70 /im is 691.2 ± 61.1 mJy (Table 4), and is taken from Ossenkopf & Henning (1994) 
to be 56 ± 11 cm^ g^-*^ at 70 ;um, assuming a standard MRN dust distribution (Mathis et al. 
1977) and an ISM-type composition consisting of graphite and silicate grains. We derive a 
total dust mass of (9 ± 2) x 10^"^ Mq, which agrees within a factor of two compared to the 
value cited by Evans et al. (2003) of (4.8 ± 1.6) x 10~^ M©, and is approximately 2-4 times 
smaller than the dust mass predicted by equation (1). The discrepancy between the Spitzer 
and ISO calculated dust masses may be largely due to the different choice of Kj,, as the ICM 
flux densities at 70 //m agree within the errors stated (Table 4). k^, is the most uncertain 
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parameter in the dust mass estimate, as its value depends largly on composition and density 
assumptions. Therefore, we note that k^, could be up to an order of magnitude larger than 
the value we have invoked here. 

The diffuse emission from ICM dust in M15 is located approximately 17" to the west of 
the cluster core. The paucity of diffuse dust toward the cluster center, where the gravitational 
potential well is steepest, is puzzling. One possible explanation for this asymmetry is a 
collection of millisecond pulsars (PSRs) near the core of M15 (Sun et al. 2002). Seven PSRs 
are located within 17" of the core (Fig. 10), and the radiative environs associated with these 
objects could lead to destruction of dust grains by sputtering or other ablation processes 
and may also inhibit dust production in stellar winds. The PSR nearest to the ICM dust 
distribution observed in the MIPS image, PSR2129+1210F, is located on the northeast outer 
edge of the 70 //m emission (Taylor et al. 1993). If, on the other hand, each mass-losing star 
has contributed 0.15 M© to the ICM on average (assuming that each star will lose 0.3 M© 
over its entire lifetime), then we see that an ICM dust mass of 1 x 10~^ M© corresponds to 
mass lost from ~10^ stars. This suggests that the ICM is short-lived, as this many stars will 
have passed through the AGB superwind phase, defined as the phase in which the mass loss 
rate exceeds the nuclear burning mass consumption rate (van Loon et al. (1999) suggests 
M ~ 10~^ — 10~^ M© yr~^ for low-mass stars), in only pa 10^ yr, which is much shorter 
than the cluster's relaxation timescale. The ICM dust therefore cannot be expected to have 
relaxed and assumed the global shape of the gravitational potential well. If this is the case, 
the offset of the dust cloud from the center of the cluster would not come as a surprise. 

Sources of the ICM dust include the post main sequence mass-losing stars identified in 
Fig. 8. If we assume that the dust-to-gas ratio scales in proportion to metallicity during the 
superwind phase, as indicated by Marshall et al. (2004), then at the metaUicity of M15, a 
dust mass loss of 10""^'^ Mq yr~^ is expected (Becker 2000). With this mass loss rate, we 
again find that the dust has been accumulating for approximately 1x10^ years, significantly 
shorter than the time between subsequent passages of the galactic plane (Evans et al. 2003). 
This short time scale suggests that dust does not survive long in the ICM. Processes that 
could be responsible for removing dust from the ICM include ram pressure by the galactic 
halo gas, radiation-driven outfiow or photo-destruction. 

4. CONCLUSIONS 

Analysis of our Spitzer image data on the core of the globular cluster M15 show strong 
evidence for the presence of intercluster medium (ICM) dust in the cluster core, with a mass 
of (9 ±2) X 10~^ Mq and with an equilibrium temperature of ~ 70 K. This is the first secure. 
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high signal to noise detection of ICM dust in a globular cluster. Also present surrounding 
the core are populations of dusty AGB and post-AGB stars, along with the planetary nebula 
K648. Using IRS spectral data, we have observed both the [S IV] and [Ne II] fine structure 
lines in K648 and have derived abundance estimates. 

The unique capabilities of Spitzer have enabled us to identify both the interstellar dust 
and the dust producers in M15. This is surprising at such low metallicity ([Fe/H] = -2.4), 
and may have implications for dust production in the early universe. The mass of the ICM 
dust in Ml 5 suggests that it has been accummulating for ~ 10^ years, which is a factor 
of ten shorter than the time since the last galactic plane crossing. The dust mass is also 
approximately 4 times smaller than the mass predicted by Evans et al. (2003). Both of these 
results imply that such dust does not survive long compared to its production rate, and is 
thus part of a stochastic process. 

We acknowledge helpful discussions with G. Schwarz and B. Moore regarding approaches 
to abundance modeling in PNe, and BDM for providing source code for multi-level atom 
calculations. M. L. B., C, E. W., L. A. H., E. P., and R. D. G. are supported in part by NASA 
through Spitzer contracts 1256406, 1215746, 1268006, and 1276760 issued by JPL/Caltech 
to the University of Minnesota, and by National Science Foundation grant AST02-05814. 
This publication makes use of data products from the Two Micron All Sky Survey, which is 
a joint project of the University of Massachusetts and the Infrared Processing and Analysis 
Center/California Institute of Technology, funded by the National Aeronautics and Space 
Administration and the National Science Foundation. 
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Fig. 1. — Three-color images of the core of M15 (NGC 7078) spanning 3.6 /xm to 70 yum. 
The field shown is 3' x 6.3'. (a) Blue is 3.6 /xm, green is 4.5 /xm, and red is 8 /xm. The PN 
(K648) becomes visible near the upper left of the cluster core, (b) Blue is 5.8 //m, green is 
8 //m, and red is 24 //m. The intra-cluster medium (ICM) becomes visible as two resolved 
high surface brightness patches (IR2 & IRlb) off-center from the cluster core. Three dusty 
sources, IR3a, IR3b, and IR4, also become prominent at 24 /im.(c) Blue is 8 /xm, green is 
24 /im, and red is 70 /im. The two ICM clouds are no longer resolved at 70 /im (IRla). 



-23- 



CO 

E 



CO 
CO 




-0.6 -0.4 



-0.2 0.0 0.2 
K-3.6jum Color 



0.4 0.6 



Fig. 2. — Color magnitude diagram showing [K-3.6 /xm] colors. K magnitudes are from 
Cohen et al. (2005). The dashed hne corresponds to the median color, which is consistent 
with K giant stars. 
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Fig. 3. — IRAC color magnitude diagrams (CMDs) for M15. Fluxes are not color corrected 
or dereddcned. (a) CMD for 4.5 /im-3.6 /im color vs 3.6 fira flux, (b) CMD for 8 /im-3.6 fim 
color vs 3.6 /xm flux. The uncertainties in the photometry are given in Fig. 5. The regions 
enclosed by the dashed lines (left of vertical, above horizontal) are the mass loss regions. 
The locations of these stars are plotted in Fig. 8, and their fluxes are listed in Table 3. 
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Fig. 4. — MIPs color magnitude diagrams (CMDs) for M15. Fluxes are not color corrected 
or dereddened. (a) CMD for 24 ^m-3.6 fim. color vs 3.6 fira flux, (b) CMD for 24 /xm-8 /xm 
color vs 8 /xm flux (see Fig. 5 for uncertainties in the photometry). The diamonds mark 
the mass losing stars identified in Fig. 3. The dash-dot lines mark the approximate division 
between AGB and post-AGB stars (Groenewegen 2006). IR3a, IR3b, and IR4 each fall well 
within the post-AGB region in (b), while only IR3a and IR3b fall within the post-AGB 
region in (a). The IRAC colors of the other mass losing sources suggest that they are AGB 
stars. 
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Fig. 5. — Photometry errors for the IRAC 3.6, 4.5, 8, and MIPS 24 /xm bands. 
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Fig. 6. — Spitzer IRS spectrum of the PN (K648) obtained with the short low (SL) module. 
The most prominent emission line in the spectrum is [Ne II] A = 12.81 fim in addition to 
weak H recombination lines and the [S IV] line. The representative error bar [left inset) 
shows the average photometric uncertainty in the spectra. No strong dust continuum is 
evident in the IRS spectrum. 
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Fig. 7. — Three color image of 8 /xm (blue), 24 /im (green), and 70 /xm (red) images. The 
24 /im and 8 /im images were rebinned to the 70 /im image pixel size of 4.925 arcseconds 
then convolved with the 70 //m Spitzer point response function (PRF). This image shows 
that the 70 /im emission is not unresolved stellar emission. 
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Fig. 8. — Locations of Mass-losing stars in M15. (a) the 3.6 /im mosaic and (b) the 24 /xm 
mosaic. All of the stars marked are located redward of the RGB. 
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Fig. 9. — Spectral energy distribution (SED) of the stellar and dust components near the 
core of M15. Fluxes (see Table 4) were obtained by using a square aperture centered on the 
dust emission at ISO wavelengths. The best two-blackbody fit (x^ = 5.26) yields 4698 ± 
58 K and 70 ± 2 K for the stellar and dust components, respectively. The dash-dot lines are 
the individual blackbodies, while the solid hue is the sum of the two blackbodies. 
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Fig. 10. — MIPS 24 |Um mosaic overlayed with radio contours (1.4 GHz) (Condon et al. 
1998). There are ten contour levels ranging from pixel values of 1.56 to 5.0 mJy/Beam on a 
square root scale. Pulsar locations are marked with an "+" symbol, and the center of the 
cluster is marked by a diamond. The inset (dotted lines) has dimensions of approximately 
1.3' X 1.0'. 
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Table 1. M15 properties^ 



[Fe/H] 


Core radius 


Z-height 




^^cluster 




E(B-V) 


D 




(") 


(kpc) 


(years) 


(Mo) 


(km s-i) (Lo) 




(kpc) 


-2.4 


4.1 


4.79 


4x10^ 


2x10^ 


40.9 5.8x10^ 


0.11 


9.98±0.47 



^ Values reproduced from Evans et al. (2003), except E(B-V) from Schlegel et al. (1998) and distance 
from McNamara et al. (2004). 



Table 2. Observational Summary 



Instrument 


AORkey 


Exposure Time^ 
(s) 


Coverage 


pixel size 
("/pixel) 


IRAC 


12030208 


10.4 


9.5' X 9.5' 


1.22 


MIPS 24 urn 


12030464 


9.96 


11' X 12' 


1.245 


MIPS 70 iim 


12030464 


10.49 


8.5' X 14.5' 


4.925 


IRS 


15733760 


60 (SL) 







^ Exposure time of individual BCD images. For IRAC, only medium 
exposure HDR frames are listed. 



Table 3. Flux Densities^ of Dusty IR Sources 



ID'' 


r 

arcsec 


3.6 //m 


4.5 //m 


5.8 iim 
mJy 


8 iiui 




24 fj,m 


2MASS J21295231+1210515 


0, 


.212 


13.99±0.70 


10.87±0.54 


7.02±0.35 


4.22±0, 


.21 


0.44±0.02 


2MASS J21295264+1210440 


0, 


.201 


6.24±0.31 


10.21±0.51 


8.64±0.43 


5.09±0, 


.25 


0.66±0.02 


SSTU J212953.55+120910.7 


0, 


.982 


8.91±0.45 


9.28±0.46 


6.80±0.34 


3.91±0, 


.20 


0.45±0.02 


2MASS J21295383+1209338 


0, 


.200 


6.15±0.31 


1.43±0.07 


5.22±0.26 


5.00±0, 


.25 


0.61±0.02 


2MASS J21295473+1208592 


0, 


.032 


4.80±0.24 


6.87±0.34 


5.23±0.26 


2.92±0, 


.15 


0.30±0.02 


2MASS J21295618+1210179 


n 


.uou 


21.08±1.05 


24.61±1.23 


22.05±1.10 


>1.78±0, 


.09 


2.93±0.02 


2MASS J21295678+1210269 


0, 


.271 


30.65±1.53 


30.35±1.52 


>6.19±0.31 


12.86±0, 


.64 


1.25±0.02 


2MASS J21295703+1209376 


0, 


.242 


4.79±0.24 


23.99±1.20 


17.84±0.89 


9.78±0, 


.49 




2MASS J21295712+1210043 


0, 


.190 


11.11±0.56 


1.73±0.09 


4.79±0.24 


5.08±0, 


.25 




2MASS J21295716+1209175 


0, 


.250 


17.59±0.88 


15.18±0.76 


12.41±0.62 


3.56±0, 


.18 


0.67±0.02 


2MASS J21295758+1209552 


0, 


.301 


12.34±0.62 


12.00±0.60 


>0.27±0.01 


3.45±0, 


.17 





Table 3 — Continued 



ID'' 


r 

arcsec 


3.6 /im 


4.5 /im 


5.8 /im 
mJy 


8 /im 


24 iim 


2MASS J21295756+1210276 


0.022 


6.06±0.30 




>0.62±0.03 


2.55±0.13 


0.44±0. 


,02 


2MASS J21295815+1209466 


0.493 


18.20±0.91 


28.30±1.41 


21.12±1.06 


8.83±0.44 


1.66±0. 


,02 


2MASS J21295828+1209280 


0.202 


4.85±0.24 


4.13±0.21 


0.91±0.05 


0.33±0.02 




- 


2MASS J21295832+1209128 


0.181 


8.03±0.40 


7.32±0.37 


5.40±0.27 


4.98±0.25 


0.50±0. 


,02 


2MASS J21295881+1209285 


0.301 


9.78±0.49 


10.12±0.51 


5.39±0.27 


2.27±0.11 


0.47±0. 


,02 


2MASS J21295937+1210029 


0.420 


8.35±0.42 


7.35±0.37 


8.61±0.43 


1.06±0.05 






SSTU J212959.69+120739.0^ 


4.722 


0.31±0.04 


0.29±0.02 


0.26±0.02 


1.47±0.07 


2.37±0. 


,02 


2MASS J21295981+1211107 


0.231 


6.41±0.32 


5.73±0.29 


4.14±0.21 




0.28±0. 


,02 


2MASS J21295996+1207282'^ 


0.100 


0.34±0.03 


0.25±0.02 


0.24±0.02 


2.05±0.10 


2.52±0. 


,03 


2MASS J21300062+1209284 


0.140 


2.95±0.15 


4.05±0.20 


>0.12±0.01 


1.94±0.10 


0.22±0. 


,02 


2MASS J21300097+1210375 


0.054 


5.78±0.29 


6.11±0.31 


4.41±0.22 


2.50±0.12 


0.31±0. 


,02 



Table 3 — Continued 



ID^ 


r 


3.6 /im 


4.5 /im 


5.8 /im 


8 /im 


24 //m 




arcsec 






mJy 






2MASS J21300277+1206557^ 


0.010 


9.34±0.47 


6.18±0.31 


4.04±0.20 


3.29±0.16 


3.66±0.02 



'^Flux densities are color corrected and dereddened using E(B-V) = 0.11 and A;^ values from Indebetouw et al. 
(2005). 

'^Standard 2MASS designations arc listed for all stars that have a 2MASS counterpart located within r = 0.5" from 
the Spitzer coordinates. Otherwise, stars are listed with standard Spitzer designations where U stands for "uniden- 
tified". 

"^Designated as lR3a in text. 3.6, 4.5, and 5.8 /imfluxes were determined with aperture photometry in IDL, as 
this source was not detected by APEX at these wavelengths. 

"^Designated as IR3b in text. 3.6, 4.5, and 5.8 //mfluxes were determined with aperture photometry in IDL, as 
this source was not detected by APEX at these wavelengths. 

^Designated as IR4 in text. 
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Table 4. ICM Dust (IRla) Flux Densities 



A 


Fa" 


Instrument (s) /Mission(s) 




mJy 




0.56 


4667.±233. 


Hbl Wl'rC2° 


1.25 


10857.±186. 


2MASS 


1.65 


10450.±204 


2MASS 


2.16 


7342.±134. 


2MASS 


3.6 


3669.±1. 


IRAC 


4.5 


2336.±1. 


IRAC 


5.8 


1442. ±1. 


IRAC 


8.0 


757.7±0.5 


IRAC 


24 


159.4±0.1 


MIPS 


60 


515.±155. 




70 


691.2±61.1 


MIPS 


70 


578.±173. 




90 


303.±91. 





^ Aperture size is 130.5 x 130.5 arcseconds, cor- 
responding to the ISO aperture from Evans et 
al. (2003). All fluxes are dereddened with E(B- 
V)=0.11. K\ values come from Rieke & Lebofsky 
(1985) and Indebetouw et al. (2005). 

^ From image U5B80101R from the pubhc HST 
archives. 



From Evans et al. (2003). 
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Table 5. Spitzer Infrared Line Fluxes of K648 



A 








/im 


Ion/Line 


ergs cm~^ s~^ 




12.81 


[Ne II] 


(1.74 ± 0.08) X 10- 


-13 


12.35 


HI 7-6 


(9.75 ±3.08) X 10" 


-15 


10.51 


[S IV] 


(8.24 ±3.56) X 10" 


-15 


7.46 


HI 6-5 


(4.38 ± 0.68) X 10" 


-14 



^ Integrated line fluxes assuming a fixed IRS 
instrumental resolution of 0.0404 /xm, and a 
Gaussian line profile. 
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Table 6. Sulfur and Neon Abundances for K648 



Quanity 


Value 




Te(K) 




iz,oUU 


nel^cm ) 




i, ( UU 


jQb IIJo717,D7dlj/i(Hpj 


(i.by ± i.izj 


X 10 


Q+ /"tJ+ 


( on _i_ o nnN 
(o.c)U ± o.UUj 


X lU 


J([o llljyoozj/J (rlpj 


/'y c^n _i_ o otN 


1 n — 3 
X lU 


Q2+ /■tT+ 


(z.oo ± i.ioj 


1 n— 8 
X iU 


iVJiU.oi /im/j(rlpj 


[pAl ± z.o4j 


X iU 


C3+ /XT+ 

b / il^ 


/I in _L n o\ 
(i.iU ± U.ozj 


X lU 


iotal o/H 


/ A OO _L 1 00\ 

(^4.zo ± i./o j 


1 n— 8 
X iU 


/([Ne III]3967)//(H/?)'i 


(2.7 ±0.9) 


X 10-3 


Ne2+/H+ 


(3.69 ± 1.23) 


X 10-6 


7([Ne IIl]3969)//(H/3)'i 


(1.19 ±0.30) 


X 10-^ 


Ne2+/H+ 


(4.89 ± 1.23) 


X 10-6 


7([Ne II] 12.81 //m/7(H/3)^ 


(1.15 ±0.18) 


X 10-^ 


Ne+/H+ 


(1.53 ±0.21) 


X 10-5 


Total Ne/H 


(2.39 ±0.27) 


X 10-5 




(5.0 ± 1.3) 


X 10-5 


[S/0]f 




-2.64 


[Ne/0] ^ 




±0.54 



* From Barker ct al. (1983). 

^ Prom Barker et al. (1983). Note lines are blended. 

H/3 flux derived from Spitzer observations of H I re- 
combination lines (see §3.2.1). 

Prom Adams et al. (1984). Measured values corrected 
for extinction (see §3.2.2). 

^ Prom Pena et al. (1993). 
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logw{Xobject) - logio{XQ). 



This figure "fl.png" is available in "png" format from: 



http://arXiv.org/ps/astro-ph/0606236vl 



This figure "fS.png" is available in "png" format from: 



http://arXiv.org/ps/astro-ph/0606236vl 



This figure "fT.png" is available in "png" format from: 



http://arXiv.org/ps/astro-ph/0606236vl 



This figure "fS.png" is available in "png" format from: 



http://arXiv.org/ps/astro-ph/0606236vl 



